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1 .  INTRODUCTION 


The  recent  development  of  various  kinds  of  tunable  solid-state 
lasers1”3*  has  revitalized  interest  in  electronic-vibrational  (vibronic) 
transitions  of  impurity  ions  in  solids;  in  particular,  laser  action  on 
phonon-terminated  transitions,  predicted  by  McCumber,H  has  been  real¬ 
ized.  A  thorough  understanding  of  the  mechanism  of  these  phonon- 
terminated  lasers  can  and  should  lead  to  the  development  of  a  useful  and 
versatile  class  of  lasers  that  operate  at  shorter  and  at  longer  wave¬ 
lengths  than  do  presently  existing  systems. 

The  mechanism  for  vibrationally  induced  electric-dipole  transitions 
for  transition-metal  ions  at  a  site  with  inversion  symmetry  in  a  solid 
was  discussed  by  Van  Vleck3  as  early  as  1937;  however,  serious  investi¬ 
gations  into  the  order  of  magnitude  of  vibronic  electric-dipole  transi¬ 
tions  did  not  begin  until  the  1950's.  During  this  latter  period,  a 
large  number  of  papers  were  written  invoking  vibronics  in  the  analysis 
of  transition-metal  ions  with  an  unfilled  d  shell.  Of  these  papers,  the 
work  of  Liehr  and  Ballhausen^ ' 3  purported  to  be  the  most  complete, 
although  from  our  standpoint  a  further  discussion  of  details  would  have 
been  helpful.  An  early  paper  by  Koide  and  Pryce9  has  been,  by  far,  the 
most  useful  to  us  in  our  efforts. 

Simultaneously  with  intense  research  efforts  on  the  role  of  vibron¬ 
ics  in  the  behavior  of  transition-metal  ions,  some  work  was  also  begun 
on  the  rare-earth  series.  In  two  papers,9”113  Satten  discussed  the 
effects  of  vibronics  on  c lectric-dipole  transitions  and  the  Zeeman 
splittings  of  rare-earth  ions  in  solids  and  solutions.  Also,  Pollock 
and  Satten11  and  Satten  et  al12  extended  the  analysis  to  actinide  ions 
as  well  as  rare-earth  ions.  Much  of  the  latter  work  relied  very  heavily 
on  group-theoretical  methods  in  the  development,  in  order  to  extract  the 
symmetry  dependence  very  concisely.  Much  of  this  analysis  was  applied 
to  cases  where  the  ions  under  investigation  occupied  sites  of  low  symme¬ 
try  or  cubic  sites  in  solids  where  charge  compensation  is  required  (for 
example,  triply  ionized  rare-earth  ions  in  CaFj). 

The  discovery13  of  the  elpasolite  hexachlorides,  Cs2NaRClg  with  R 
any  of  the  triply  ionized  rare-earth  ions,  changed  the  experimental 
situation.  It  was  found  that  CS2NaRClg  of  good  optical  quality  could  be 
grown  with  R  any  rare  earth,  or  with  R  =  Y  and  a  small  amount  of  rare 
earth  added.  The  site  occupied  by  R3+  in  CsnNaRClg  has  0^  symmetry,  and 
the  resulting  spectra  are  rich  in  vibronics.14  Because  of  the  sharpness 
of  the  rare-earth  spectra,  a  large  number  of  different  vibrational  modes 
can  be  identified,  and  by  fitting  the  frequencies  of  these  modes  a 


’•'Because  of  their  large  number,  literature  references  are  listed  at 
the  end  of  the  text  rather  than  on  each  page. 
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number  of  force  constants  in  the  vibrational  Hamiltonian  can  be  deter¬ 
mined.  18  The  determination  of  these  force  constants  allows  a  quantita¬ 
tive  calculation  of  the  intensity  of  the  vibrationally  induced  electric- 
dipole  transitions.  Such  a  calculation  has  been  performed  by  Hasan  and 
Richardson,18  whose  theoretical  predictions  were  found  to  be  in  rough 
agreement  with  experiment.  These  calculations  are  important  for  possi¬ 
ble  technical  applications17  of  this  material  or  the  less  hygroscopic, 
isomorphic  compounds  CSjNaRFg  and  RbjNaRFg.  Unfortunately,  very  few  of 
the  vibronic  spectra  of  the  latter,  more  practical  compounds  have  been 
reported.  Nevertheless,  an  analysis  of  vibronics  of  Er3  +  in  Cs2NaRClg 
should  be  very  useful. 

In  this  report,  we  investigate  the  electric-dipole  transitions  of 
Er3  +  in  Cs0NaErClg  caused  by  the  odd  vibrations  of  the  (RClg)-3  complex. 
Most  of  the  results  will  be  presented  in  a  form  that  allows  simple 
generalization  to  the  complete  unit  cell  (in  the  k  =  0  limit),  although 
this  generalization  will  necessarily  involve  the  use  of  a  computer. 
Thus,  the  analysis  of  the  normal  modes  of  the  complete  unit  cell  such  as 
given  by  Lentz18  can  be  readily  used  in  any  future  calculations.  In 
section  2,  we  discuss  the  modes  of  vibration  of  the  (RClg)  8  ion.  In 
section  3,  the  static  crystal  field  is  discussed  using  the  point-charge, 
point-multipole  potential  tTiodel.  Also  in  section  3,  the  self-induced 
crystal  field  of  Faulkner  and  Richardson18  is  presented  in  a  slightly 
different  form  generalized  from  dipoles  to  multipoles.20  The  vibration¬ 
al  electronic  interaction  is  the  subject  of  section  4,  along  with  an 

explicit  calculation  of  the  electric-dipole  transitions  for  Er3+  accom- 

_  "} 

panying  the  Vg  vibration  of  (ErClg)  . 

2.  VIBRATIONAL  MODES  OF  (RClg)-3 

The  ion  cores  that  make  up  the  crystal  and  determine  the  crystal 
field  felt  by  an  electron  on  a  rare-earth  ion  are  located  at  sites  in 
the  crystal  which  are  local  minima  of  the  effective  ion-core  potential; 
i.e.,  each  ion  is  in  its  own  potential  well.  It  is  well  known  that 
these  ions  can  execute  small  vibrations  within  their  wells,  and  that  for 
small  amplitudes  of  vibration  each  ion  behaves  like  a  harmonic  oscilla¬ 
tor.  The  interionic  forces  will  in  general  couple  these  oscillators,  so 
that  the  crystal  as  a  whole  has  extended  wave-like  normal  modes  which, 
for  an  infinite  crystal,  can  be  characterized  by  a  wave  number,  k.  When 
these  modes  are  quantized,  their  quanta  make  up  the  elementary  excita¬ 
tions  of  the  background  crystal — the  so-called  phonons. 

The  normal  nodes  of  a  typical  crystal  fall  into  two  categories: 

(1)  Acoustic  modes:  these  are  characterized  by  all  the  atoms  of  a 
unit  cell  moving  at  once  in  nearly  the  same  direction  by  nearly  the  same 
amount.  Every  crystal  has  three  such  modes. 
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(2)  Optical  modes:  these  modes  involve  some  countermotion  of  atoms 
in  the  unit  cell;  if  there  are  N  atoms  in  the  unit  cell,  there  are  N  -  3 
optical  modes. 

In  general,  acoustic  modes  correspond  at  low  values  of  k  to  crystal 
strain  waves,  and  involve  macroscopic  motion  of  the  crystal  bulk;  their 
frequencies  u>(k)  cover  a  rather  broad  band  in  k-space:  from  k  =  0  to  k 
on  the  order  of  a  reciprocal  lattice  vector,  u(k)  may  vary  by  as  much  as 
100  meV.  By  contrast,  the  optical  modes  are  always  microscopic  motions 
which  do  not  give  rise  to  bulk  movement  of  the  crystal  matrix;  their 
frequencies  tend  to  be  narrow  band  (~5  meV)  due  to  the  rather  weak 
coupling  of  adjacent  unit  cells  that  execute  this  kind  of  motion.21 
Because  the  form  of  the  crystal  field  is  intimately  related  to  the  unit 
cell  configuration,  the  motions  associated  with  these  two  kinds  of 
normal  modes  will  have  radically  different  effects  on  the  electronic 
states  of  a  rare-earth  ion. 

In  dealing  with  a  crystal  with  a  complicated  unit  cell  like  elpaso- 
lite,  it  pays  to  examine  the  crystal  vibrations  from  a  molecular  stand¬ 
point.  This  entails  the  following  simplified  picture  of  the  ion  motion: 

(1)  Acoustic  modes  correspond  to  simple  translation  of  every  ion  in 
the  cell,  and  have  zero  frequency. 

(2)  Optical  modes  do  not  couple  one  unit  cell  to  another. 

Clearly  this  picture  is  simply  one  of  a  single  unit  cell  hanging  in 
space,  and  being  treated  like  an  isolated  molecule.  We  expect  that  this 
picture  will  be  fairly  good  for  those  normal  modes  which  are  self- 
contained,  i.e.,  those  for  which  the  boundary  atoms  of  the  cell  are 
motionless.  Actually,  it  turns  out  to  be  good  for  all  the  optica? 
modes,  with  three  glaring  exceptions:  an  isolated  unit  cell  has  three 
rotational  modes  with  zero  frequency  which  go  into  nonzero  frequency 
optical  modes  that  will  clearly  couple  the  unit  cells  together  in  a 
nontrivial  way. 

In  figure  1  we  show  the  full  unit  cell  of  elpasolite.  In  general,  a 
full  normal-mode  analysis  of  the  optical  modes  of  this  crystal  reveals 
that  one  need  only  specify  the  motion  of  the  atom  in  the  primitive  unit 
cell,  and  that  translation  operators  will  give  the  motion  of  the  remain¬ 
ing  atoms.  These  essential  atoms  (from  a  vibrational  standpoint)  are 
shown  in  figure  2,  taken  from  a  paper  by  Lentz.18  The  basic  vibrating 
structure  consists  of  a  pair  of  octahedra  joined  at  a  vertex,  flanked  by 
two  cesium  atoms.  It  is  clear  that  if  a  rare  earth  is  placed  at  the 
center  of  one  of  these  octahedra,  the  normal  modes  which  couple  most 
strongly  to  it  are  those  which  correspond  to  the  vibrations  of  the  free 
octahedral  complex. 
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Figure  2.  Vibrational  unit  cell  (after  ref  18;  equivalent  atoms 


have  the  same  labels). 
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The  normal  modes  of  a  free  octahedron  have  been  known  for  some 
time.®  In  figure  2  we  show  an  octahedron;  labeling  the  atomic  dis¬ 
placement  coordinates  (xN,  yN,  zN),  N  =  1 ,  .  .  .  7 ,  we  see  that  there 
are  21  degrees  of  freedom,  and  hence  21  normal  modes.  A  group-theoretic 
analysis  reveals  that  these  21  modes,  which  form  a  representation  of  the 
0^  group,  can  be  grouped  by  symmetry  into  the  following  irreducible 
representations  of  0^: 


I\(a,  ) 

i  i  g 

r!(e,  ) 

3  ig 

r*(t,  ) 

4  1  g 

3IYSu’ 

r5<t2U> 


W 


Of  these  representations,  the  t1g  and  one  of  the  tlu  modes  correspond  to 
the  zero-frequency  rotation  and  translation  modes,  respectively. 


Now,  there  exist  sets  of  so-called  "symmetry  coordinates,"  i.e., 
linear  combinations  of  the  (xN,  yN,  zN),  which  reduce  the  classical 
vibrational  Hamiltonian  to  block  form:  by  this  we  mean  that  the  symme¬ 
try  coordinates  transform  as  irreducible  representations  of  0^.  If  a 
given  irreducible  representation  appears  once  in  the  regular  representa¬ 
tion  generated  by  the  ion  coordinates  (xN,  YN,  zN),  then  the  correspond¬ 
ing  symmetry  coordinate  is  a  normal  mode.  In  general,  however,  the 
regular  representation  may  contain  a  given  irreducible  representation 
several  times;  in  this  case,  there  will  be  several  symmetry  coordinates 
associated  with  the  ion  coordinates,  none  of  which  is  a  normal  mode.  In 
this  case,  the  symmetry  coordinates  do  not  fully  diagonalize  the  Hamil¬ 
tonian,  and  the  true  normal  coordinates  are  linear  combinations  of 
symmetry  coordinates  of  a  given  irreducible  representation. 

The  problem  of  a  free  vibratinq  octahedral  complex  illustrates 
everything  we  have  discussed  thus  far.  Figure  3  shows  the  lahelinq 
system  that  we  will  use  here.  Fiqure  4  shows  the  eiqht  "svmmetrv  vibra¬ 
tions"  of  an  octahedron  with  a  rare-earth  ion  at  its  center,  correspond¬ 
ing  to  the  appropriate  symmetry  coordinates;  if  the  actual  ion 
coordinates  (the  six  chlorine  ions  plus  the  rare  earth)  are  labeled  (xN, 
yN,  zN),  then  in  terms  of  these  coordinates  the  svmmetrv  coordinates  are 
as  follows: 
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(A)  Translation:  T.(t  ) 

4  1  u 

-  n2y_  4.  r  2  f  v  4.  v.  4-  V  _  4-  1 


1  3  ~ 
where 


n2z?  + 


C  =  (6  +  e2  j-1/2  ( 

e  =  /Mq/M  , 

H  "Cl,  t 

Mq  is  the  rare-earth  mass,  and 
M  is  the  ligand  (i.e.,  chlorine)  mass. 


Fiqure  3.  Labeling  for 
RClg  octahedron. 


Note  that  the  are  simply  the  center-of-mass  coordinates  of  the  octa¬ 
hedral  complex. 


(B)  Rotation:  T.(t  ] 
4  1  g 
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2  ^“2 
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X2  +  yi  -  V 


(C)  "Breather:"  I\  (a,  ) 

1  1g 

4  '  *  y2 
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3  1  g 
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10 
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4  1  u 


r 

S1 

+ 
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If  we  write  the  classical  Hamiltonian  in  terms  of  these  eight  symmetry 
coordinates  (with  their  appropriate  conjugate  momenta),  it  is  completely 
diagonal  except  for  a  term 


This  term  is  not  diagonal  because  0  and  E,  both  transform  like  r^(t1u). 
The  normal  modes  depend  on  the  values  of  V^y  which  in  turn  relate  to 
the  specific  forces  between  the  octahedron's  constituents. 

If  we  now  want  to  consider  the  effect  these  vibrations  have  on  the 
crystal  field  experienced  by  the  octahedron's  center,  we  clearly  need  to 
know  the  detailed  motion  of  each  ion  when  a  given  normal  mode  is  excit¬ 
ed.  In  quantum  mechanical  language,  there  are  phonons  associated  with 
each  normal  mode;  the  motion  of  an  electron  will  cause  emission  and 
absorption  of  these  phonons,  and  these  phonons  in  turn  act  back  on  the 
electron.  The  interaction  potential  will  be  derived  in  section  4;  for 
now,  we  note  that  it  will  involve  the  ion  coordinates  themselves,  not 

the  normal  modes.  In  order  to  evaluate  the  various  transition  probabil¬ 

ities,  we  must  therefore  express  the  ion  coordinates  in  terms  of  the 
normal  modes;  i.e.,  we  must  invert  relations  (A)  to  (H).  Once  this  is 
done,  we  must  subtract  off  the  rare-earth  coordinate,  since  the  electron 
is  "riding  with"  the  rare  earth.  When  this  is  done,  the  following 
results  obtain:  if  Ax^  =  ,  Ay^  =  y^  -  y-j,  and  Az^  =  -  z 7, 

then 

Axi  =  -A61  +7?  °0  ' 

Ay)  =  -B£2  -  c®2  +  (<2  +  p3  +  v3 )  ' 

A2 1  =  _b^3  ”  c03  +  ~2^  i  -IC3  _  p2  +  v2~) 

Ax2  —  ^K-|  —  P2  +  , 

Ay2  =  -A02  +  A  -  (o,  +  02))  , 

AZ2  =  "B53  “  C63  +  2C  +  P 1  ‘  '• 
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Note  that  the  translation  mode  does  not  contribute  to  the  relative  coor¬ 
dinate.  This  can  be  traced  to  the  assumption  that  all  the  unit  cell 
atoms  move  by  exactly  the  same  amount  in  this  mode;  for  k  =  0  this  is 
true,  but  for  nonzero  k  the  acoustic  mode,  which  corresponds  to  the 
coordinate  5,  will  indeed  couple  to  the  electron's  motion,  giving  rise 
to  the  so-called  "deformation  potential." 


3.  CRYSTAL- FIELD  INTERACTION 


In  the  phenomenological  theory  of  the  crystal 
an  electron  on  a  rare-earth  ion,  the  conventional 


field  interacting  with 
Hamiltonian22 


H 


Y 

L 

nm 


B* 

nm 


C  (i  ) 
nm 


(1  ) 


is  frequently  used.  The  Bnm  are  obtained  by  fitting  the  experimentally 
measured  energy  levels  using  a  suitable  set  of  free-ion  wave  functions 
to  evaluate  the  matrix  elements  of  equation  (1).  The  number  of  Bnm  in 
equation  (1)  is  restricted  by  the  symmetry  of  the  site  occupied  by  the 
rare-earth  ion.  For  the  material  CsjNaRClg,  the  site  symmetry  of  the 
rare-earth  ion  is  O^,23  so  that  the  only  independent  parameters  are  B4f) 
and  Bg0;  the  Bn4  are  related  by 

B44  =  5B40//7B  ' 

(2) 

B64  =  -  /7/2  B60  • 

The  spectrum  of  Er3+  in  CS2NaErClg  has  been  analyzed  by  several 
workers1^ '  >25  ancj  the  resulting  Bnm  are  given  in  table  1.  As  can 
be  seen,  there  is  considerable  disparity  in  the  reported  values  of  B4q 
and  B60.  This  disparity  is  due  in  part  to  the  different  identifications 
made  of  experimental  levels,  and  in  part  to  use  of  different  free-ion 
wave  functions  in  calculations. 

The  oldest,  and  perhaps  the  most  reliable, 
crystal-field  theory  is  the  point-charge  model 
developed  by  Bethe2°  in  his  paper  on  the  spectra 
of  transition  metal  ions  in  a  solid.  In  such  a 
model,  the  constituent  ions  at  Rj  from  the  rare- 
earth  ion  nucleus  have  charges  eq^,  and  the 
rare-earth  electron  at  r  interacts  with  these 
ions  through  the  Coulomb  potential,  given  by 


TABLE  1 .  REPORTED  B4Q 
AND  B,n  (cm-1)  FOR 
ErJ  IN  Cs2NaErCl6 


B40 

03 

CTi 

O 

Ref  . 

1  384 

32 

24 

1  368 

32 

25 

1425 

153 

26 

1602 

151 

27 

15 


(3) 


„0  -  ,2  T  J 

CEF  "  l  +  + 

3  I Rj  -  H 

A  multipolar  expansion  of  equation  (3)  can  be  readily  obtained  to  qive 

HcV  -  -e2  I  I  ‘Jjrncn»(r“)C5m(«jVRj  +  1  '  “ 

j  nm 

and  if  we  sum  over  all  the  electrons,  we  may  write 


=  y  a°*  y  rnc  r;.  1 

^  nm  L  l  nm  v  1  ' 


\0m  =  -e2  I  qjCnm(^j)/Rj  +  1  *  (6) 

j 

These  A8  are  frequently  referred  to  as  the  multipolar  components  of  the 

crystal  field.  The  A8  have  been  evaluated  tor  a  large  number  of 
o  *7  nm 

solids,-^7  the  results  in  reference  27  are  given  in  such  a  way  that  the 
q^  in  equation  (6)  can  be  replaced  by  effective  charges  on  the  ligands. 
Such  a  procedure  can  be  used  to  partially  account  for  covalent  complexes 
in  the  solid. 

Hutchings  and  Ray28  recognized  that  in  solids  where  the  ligands 
occupy  sites  of  sufficiently  low  symmetry,  the  point  charges  of  an  ionic 
solid  could  induce  electric  dipoles  and  electric  quadrupoles  in  these 
ions.  These  workers  did  not  take  into  account  the  fact  that  these 
dipoles  and  quadrupoles  also  contributed  to  the  electric  fields  at  the 
sites  of  low  symmetry.  The  consistent  solution  for  dipoles  has  been 
given,28  and  recently  the  full  self-consistent  solution  through  quadru¬ 
poles  has  also  been  done.*  In  addition,  the  corresponding  self- 
consistent  solution  for  arbitrary  multipolar  moments  has  recently  been 
published . 28 

For  a  multipolar  moment  Q  (j)  at  R.,  the  A 8  of  equation  (6)  can  be 
,  .  ,  on  KB  "]  * 

generalized 


A^m  =  -e  J  (-1  )k  (22k2t  2  <n(m)k(  0)  jn+k(m+g)> 

j  >  6 


^k  6^  ^n+kjin+B^j  ) 


where 


(2K+2n  1 
1  2k  j  ~ 


[25  +  2k  )  ! 
!  2k) ! ( 2n ) ! 


*M.  Faucher  and  0 .  Malta  ( private  communication! . 
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and  the  quantity  in  angular  brackets  is  a  Clebsch-Gordan  eoeff icient . 
It  should  be  noted  that  if  k  =  0  in  equation  (7),  the  result  is  the  same 
as  that  of  equation  (6)  if  Qqq(j)  =  eq j ;  the  two  expressions  are  consis¬ 
tent.  Thus,  if  the  (j)  have  been  determined,  their  contribution  to 
the  crystal  field  is  given  by  equation  (7). 


A  further  interaction  recently  considered  is  the  self-induced  inter¬ 
action.^®  This  interaction  is  caused  by  an  electron  on  the  rare  earth 
inducing  a  multipole  moment  in  a  ligand,  which  then  interacts  back  on 
the  electron.  The  resultant  Anm  are  given  by 

.self-induced  _  ?  V  n  t- n  U  +  n ) !  f  2n~ 1  (n!  )2  .1 
Anm  -  "e  .1.  zTn!  rsnr~~  "  J 

(8) 


X 


(Rj )/Rj4+n+2  . 


The  results  of  equations  (8)  and  (7)  have  been  used  to  evaluate  the 
various  contributions  to  Anm  for  LaCl^,3®  and  it  is  shown  that  each  of 
these  interactions  contributes  significantly.  Although  we  are  here 
interested  in  the  material  CsjNaErClg,  we  give  the  results  for  LaCl^ 
from  Brockhouse^®  in  table  2  to  illustrate  the  significance  of  the 
contribution  to  the  Anm  given  by  equation  (8). 


TABLE  2.  ELECTROSTATIC  CONTRIBUTIONS  TO  Anm  FOR  LaClj 


Dipole  and  quadrupole  calculations  performed  using 
a1  =1.5  A^  and  a =  3.44  A^ 


A 

nm 

Point 

charge3 

Point 

dipole3 

Self¬ 

dipole 

Self- 

guadrupole 

A20 

-2460 

3273 

-120 

-64 

*40 

-794 

-22 

1  35 

107 

-185 

-80 

73 

77 

ReAM* 

-85 

17 

-118 

-40 

ImA64 

-63 

-39 

-99 

-34 

3from  J. 

B.  Gruber ,  R.  P. 

Leavitt,  and  C. 

A.  Morrison,  J. 

Che/n.  Phys*  74 

(1981),  3705  (ref  30). 

bIn  IaCIj  a  rotation  of  the  coordinates  about  the  c  axis  can  be  per¬ 
formed  to  make  the  Ag4  component  real.  This  rotation  must  be  performed 
after  the  various  contributions  to  Am  are  accounted  for. 
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The  results  given  in  equations  (6),  (7),  and  (8)  represent  a  number 
of  interactions  of  electrons  with  a  rare-earth  ion  in  a  solid  from  the 
point-charge  or  point-multipole  viewpoint.  A  number  of  other  possible 
interactions  have  been  ignored,  such  as  exchange,  covalency,  configura¬ 
tion  interaction,  etc.^1  In  our  discussion  of  the  interactions  with  the 
lattice  vibrations,  we  shall  restrict  ourselves  to  the  interactions 
given  in  equations  (6),  (7),  and  (8).  For  an  excellent  evaluation  of 
the  different  interactions  as  of  1971,  the  reader  should  consult 
Newman's  article.31 


4.  THE  VIBRATIONAL  INTERACTION  POTENTIAL 

The  interaction  potential  we  shall  consider  here  is  that  derived  by 
considering  small  displacements  from  the  equilibrium  configuration  of 
the  nucleus  and  ligands.  The  particular  potential  we  shall  consider  in 
detail  is  the  point  charge  and  its  extension  to  the  point-multipole 
potential.  We  shall  further  assume  that  the  various  charge  distribu¬ 
tions  are  not  distorted  by  their  small  excursions  from  equilibrium, 
i.e.,  the  so-called  "rigid-ion"  model.  An  alternative  interaction,  the 
superposition  model,31  has  recently  been  considered  by  Kennedy32  in  his 
analysis  of  the  effect  of  an  external  strain  on  spectra  of  rare-earth 
ions  in  CaF2< 

Frequently  the  interaction  of  an  electron  on  a  rare-earth  ion  with 
its  surroundings  is  derived  from  the  expansion33 


where  the  are  the  normal  modes  of  the  lattice.  The  3V/9Q^  in  equa¬ 
tion  (9)  represent  the  interaction  of  the  electron  with  the  normal  mode 
Q^.  In  general,  equation  (9)  is  difficult  to  use  in  detailed  numerical 
computation  but  is  excellent  for  theoretical  purposes  when  one  is  inter¬ 
ested  only  in  the  analytical  consequences  of  the  effect  of  vibrations.33 
Since  we  are  interested  in  numerical  values  of  the  strength  of  the 
electron-lattice  interaction,  we  shall  use  a  somewhat  different  expan¬ 
sion,  which  we  present  in  the  following. 

Consider  a  rare-earth  ion  at  the  origin,  an  electron  at  r,  and  a 
ligand  at  ^  that  is  displaced  by  I j .  The  ligand  is  assumed  to  have  a 
set  of  multipole  moments  Q^q*  The  electric  potential  at  the  electron  is 
given  by 
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1 


t  -r  -r  ~r  j  , 

where  R  =  K  -  r  +  s.  .  Usinu  the  two-center  expansion-'  ot 

1  3  3  ’ 

K+1 


we  obtain 

ckqK  j 


k  +  1-  -  \  (-1)d  ^  cx)k  (ej  )  |  d+K  (  ota-t^  )  >  s 


Ckq[R1 )/R1  ' 


a  a 


c 

x - 


aa^sj  ^a+k  ,q  +  a ^ R2  ) 


R1 


a+k  +  1 


(ID 


where  =  R_.  -  r .  A  repeat  ot  the  two-center  expansion  ot 


Ca%k,q  +  >2VR! 


a+k  +  1 


gives 


[r2) 

_ l_~ _ _  V  f2n+2k  +  2au/2 

.k+a  +  1  t  ^  2n  ' 


La +k , q+a ( ^2 


x  rnC  (r  1 
nm  i  ‘ 


<n  (m  )k+a  (q  +  a)  |  n+k+a  (m+q  +  a)  > 
L"n+k+a,m+q  +  a^Rj  ) 

7,TiTk+aTr~ 

K3 


(12) 


We  are  interested  only  in  terms  linear  m  the  displacement  (a  =  1  in 
equation  (11)),  so  that  combininq  equations  (12)  and  (11)  into  equation 
(10)  qives 

*  =  -  l  (-nk  (2k+2)^2  <1  (a)k(q)|k+l  (q  +  a)>  (2n+2k+2)i/2 


kq 

a,  n,  m 


x  <n (m )k+1 (q+a) | n+k+ 1 (m+q+a)  >  rnCnm(r) 


(13) 


(1  )s jci a(sj  )Cn+k+1  ,  m+q+a ^R;j  ^ 

x  “n+k+2 

R1 

The  result  given  in  equation  (13)  is  the  potential  for  the  ligand  at 
Rj  ;  the  vibrational  interaction,  V',  in  the  Hamiltonian  is  given  by 
-eiji,  which  we  write 

V'  =  J  A*  rnc  (r)  (14) 

^  nm  nm 
nm 

by  analogy  to  the  interaction  of  the  electrons  on  a  rare-earth  ion  in 
the  absence  of  lattice  vibrations.  The  prime  on  V  in  equation  (14) 
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indicates  that  only  linear  displacements  are  considered,  rf  we  had 
considered  the  term  a  =  0  in  equation  (11),  we  would  have  had  the  inter¬ 
action  of  the  electron  in  the  absence  of  lattice  displacement,  that  is, 
the  crystal-field  interaction  given  by  equation  (7).  Summing  equation 
(13)  over  all  the  ligands  and  comparing  with  equation  (14)  gives 


nm 


k,g,a 


(-1)k^2V/2  ( 


\  2/ 


2n+2k+2 
2n  / 


1  /  2 


<  1  ( a)k (q  )  |  k+ 1  (q  +  a)  > 


x  <n (m )k+1 (q+a) | n+k+1 (m+q+a) >  (15) 

sjC*a(sj  )Cn+k+1,m+q  +  a^j  ) 

X  ®kq*3^  Rn+k+2 

3  3 

The  result  given  in  equation  (15)  can  be  cast  into  a  more  useful  form  by 
using  explicit  expressions  for  the  Clebsch-Gordan  coefficients  to 
obtain3  5 

(2k 22)1/2  <1 (a)k(q) |k+1 (q+a)>  f2n+2n+2)1 <n (m )k+1 ( q+a) | n+k+1 (m+q+a) > 

(n+k+1 +m+q+a ) ! (n+k+1 -m-q-a ) ! 

(n-m) !  (n+m  ) !  ( 1  -a) !  ( 1  +  a)l  (k-q) !  {k+q ) ! 

Then,  from  equation  (15),  we  have 

k  f  (n+k  +  1 +m+q+a  )!  (n+k  +  1 -m-q-a  ) ! 

Anm  “  e  L(n-m) !  (n+m) !  { 1-a) !  ( 1+a) !  (k-q) !  (k+q) ! 

x  l  Q*q(j)  , 

which  is  the  form  of  the  interaction  that  we  find  the  most  convenient  to 
use.  For  a  particular  normal  mode,  the  have  been  specified  in  sec¬ 
tion  2,  and  the  Q^q(j)  will  also  be  considered  known.  Thus  the  sum  over 
j  given  in  equation  (17)  can  be  performed  for  any  vibrational  mode  of 
interest,  and  the  resulting  A  for  that  particular  mode  can  be  deter¬ 
mined.  Thus,  the  result  given  in  equation  (17)  is  particularly  adapt¬ 
able  for  calculations  using  a  computer.  The  (j)  can  be  found  from  an 
appropriate  lattice  sum  and  some  self-consistent  method.  The  dipole 
terms  (Q-|q)  have  been  obtained  for  a  number  of  solids,25'36'37*  and  in 
particular  for  the  material  Cs2NaRClg25  for  R  any  of  the  rare-earth 
ions.  The  itj  in  equation  (17)  are  obtained  from  the  x-ray  data,  so  that 
with  these  data  the  sum  given  in  equation  (17)  can  be  coded  for  computa¬ 
tion  to  cover  as  many  unit  cells  as  is  necessary  for  convergence.  In 
computing  the  sums,  the  quantities  sjcia(sj)»  ^kq^^'  anc*  need  only 
be  specified  for  a  unit  cell,  and  lattice  displacements  can  then  be  used 

*M.  Fauchec  and  0.  Malta  (private  communication) . 
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(17) 


1/2 


(16) 
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to  cover  the  complete  solid  (in  the  K  =  0  limit).  It  is  also  possible 

to  approximately  „take  into  account  some  of  the  phonon  dispersion  by 

replacing  siC'1g^si^  ^  s^C)a(sj)  cos  and  evaluating  the  sum  on  j 

in  equation  (17)  for  particular  values  of  K.  However,  in  the  following 

we  shall  specialize  equation  (17)  to  the  particular  (RClg)-3  complex  and 

obtain  specific  algebraic  expressions  for  A  for  each  of  the  modes. 

nm 


5.  THE  ELECTRON-VIBRATION  INTERACTION  FOR  (RClg)  3 

_  O 

As  shown  in  section  2,  the  vibrational  modes  of  the  (RCl^)  complex 
decompose  into  the  following  representations  of  the  0^  group  (neglecting 
the  translation  and  rotation  modes) 


r?<a,  ) 
1  1  g 


2I\(t  ) 

4  1  u 


+  IVt2u) 


r!(e  ) 
3  g 


r!(t  > 

5  2g 


Here,  we  are  only  interested  in  the  vibrationally  induced  electric- 

dipole  transitions,  which  are  of  odd  parity  (T  ).  Since  the  total 

interaction  Hamiltonian  must  be  of  even  parity,  we  then  require  only  the 

odd-parity  vibronic  modes.  That  is,  we  shall  discuss  only  the  two 

modes  r.(t,  )  and  r,(t„  ).  The  latter  will  be  discussed  first. 

4  1 u  5  2u 

5.1  A  for  IV  (t  )  Mode 
nm  5  2u 


The  appropriate  vibration  for  the  single  r,-  mode  is  given  in 
figure  4(H).  The  displacements  corresponding  to  the  coordinate  are 


sjC1«(sj)  =  -216a,0  ' 

for 

j 

=  1  and 

6  , 

sjC1a^j)  =  z16a,0  , 

for 

j 

=  2  and 

4  , 

where  z1  is  the  displacement  of  the  ligand  labeled  1  in  figure  2;  these 
displacements  follow  from  the  mode  equations  given  in  section  2. 


The  position  and  multipole  moments  can  all  be  related  to  the 
ligand  at  1  by  rotations  about  the  z  axis.  That  is, 

e{,q(2'Cn+kf1,4"2)  =  eimV^,q<1)Cn+k+1,mtq(«l) 

with  a  similar  expression  for  the  relation  of  6  to  1,  which  we  denote 
symbolically  as 

(6)  =  e 11,111  ( 1  ) 


and 


(4) 


-imn/2 

e 


1 1  ) 


(19) 
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Using  equations  (18)  and  (19)  in  the  sum  over  a  and  j  in  equation  (17), 
we  obtain 


L  SjCTa^j)^q^)Cn+k.1,m+q(^) 

'■*  J 

'  -4z,ekq(',C„t|t„..-K,(N)5.,4p*2 

wi th  p  =  0 ,  tl ,  ±2  .  .  . 


(20) 


In  obtaining  the  result  given  in  equation  (20),  we  have  used 
the  fact  that  all  1 1? ^  j  have  the  same  value  (the  x-ray  determined  dis¬ 
tance  from  the  rare-earth  ion  to  the  Cl”).  Using  the  results  of  equa¬ 
tion  (20)  in  equation  (17)  for  the  monopole_  term  (for  k  =  0,  Qqq(J)  = 

eg  i),  we  obtain  the  following  nonvanishing  A 
*cl  nm 


with 


A32(0)  =  K3,/12  c42(1)  , 

A52(0)  =  K5/l2  c62(n  - 

A72(0)  =  K?/60  C  (1 )  , 

A76(0)  =  K?/i8  c86(1)  . 

Kn  =  -^clV1*"*2  ' 


(21  ) 


(Kn  will  be  used  for  the  multiplication  factors;  in  general,  its  value 
in  each  particular  case  will  be  given.)  Also,  we  have  taken  equation 
(17)  as 


A  =  I  A  (k)  , 

nm  ,  nm 
k 

where  k  represents  the  multipole  moment  on  the  ligand. 


(17a) 


For  the  dipole  terms  (Q^fj))  the  site  occupied  by  the  chlorine 
ion  (site  1  )  has  only  the  components  with  q  =  ±1  ,-  using  this  and  the 
result  given  in  equation  (20)  in  equation  (17)  gives 


A32(1)  =  K3  /24  [/?4  Q*,  C53(  1  )  -  /3  Qn  C51(1)]  , 

A52(1)  =  K5  / 480  [/I  Q*t  C73(1)  -  Qn  C71(1)]  , 

A72(1)  =  K?  /T20  [/33  Q*,  C93(1)  -  /T4  Qn  C91(1)]  , 

V1*  -  K7  /Ti?  V™  of,  S3(1)  -  211  C91(1)l  ' 


(22) 
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1 


with 


K  =  +4e2z  /Rn+3 
n  1 ' 


and  z-|  is  the  same  as  in  equation  (21) 


5.2  A  for  r  <t  )  Modes 
nm  4  1  u 


The  appropriate  normal  modes  for  T  (t  )  vibration  interaction 

4  1  u  ■*  f 

are  given  by  linear  combinations  of  the  two  vibrations  labeled  £  and  0, 

shown  in  figure  4  (E)  and  (F).  Thus,  it  is  appropriate  here  to  discuss 
the  fields  associated  with  the  two  corresponding  r  (t  )  symmetry  coor¬ 
dinates.  Later  we  will  discuss  the  selection  of  the  appropriate  combi¬ 
nations  of  these  two  symmetric  functions  that  form  true  normal  modes  of 
vibration.  Also,  since  the  Anm  for  the  other  two  components  of  these 
two  modes  are  trivially  related,  we  shall  give  only  the  coefficients  for 
the  vibrations  depicted  in  figure  4  (E)  and  (E). 

From  the  figures,  the  displacements,  s • ,  of  the  ligands  corre¬ 
sponding  to  £3  are 


sjC1cx^j)  =  215 


ctO 


and 


sjC1a(sD)  =  z26aO 


for  j  =  1 ,  2,6,  and  4 


for  3=3  and  5 


(  23) 


(  24) 


For  the  93  mode,  we  wil1  replace  z,  and  z2  by  z^  and  z^.  The  sum  on  a 
and  j  in  equation  (17)  becomes 


sjC1o/sj  )Cn+1  ,m+a^j  . 


oil 


(25) 


4z1Cn  +  1  ,m(^l)6m,4p  +  2z2Cn+1  ,m^3^  6n,  2k+l  6m,  0  ' 


where  p  =  0,  ±1 ,  ±2 


.  and  k  =  0,  1 ,  2 


Using  equation  (25)  in 


equation  (23)  results  in  the  following  A  (0)  for  the  two  modes: 

nm 


A10(0)  =  2M2C20(1)  +  Z2C20(3>: 

(0)  =  2KJ2Can(1)  +  27cun(3)i 


30' 


A50(0)  =  3K5[2C60(1)  +  Z;?Cfi0(3)] 


*2  60' 


A54(0)  =  2/5  Ks[cfi4(1)]  , 


(26) 
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A70(0)  =  4K7[2C80(1)  +  Z2C80(3) 
A74(0)  =  4/3  K7[c84(1 )]  , 
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L 


where 


Z2  =  Z2/Z1  ' 

„  2  .  n+2 

Kn  =  4e  qcl2,/R 


The  effective  A  (1)  for  the  T. (t.,,)  modes  can  be  obtained  in  a 
similar  manner  as  the  Anm(0)  for  this  mode.  Because  of  the  symmetry  or 
the  Cl"  site  in  Cs2NaRCl6  ( C4 v ) ,  the  dipoles  p1  ( j )  located  at  sites  1, 
2,  4,  and  6  are  nonvanishinq  for  a  =  ±1  only,  while  the  dipoles  at  3  and 
5  have  a  =  0  only.  From  equation  (17)--see  also  equation  (I7a)--we  have 


Anm(1) 


_p  v  (■  (n+2+m+q ) ! ( n+2-m-q ) i  1  / 2 

Rn  +  3  ~  ^  ( n-m )  !  ( n+m  )  !  (  1  -q  )  !  (  1  +q  )Tj 

x  1  ^1q(3  )sjC1  0^sj  )Cn  +  2,m+q1-Rj  ) 

j 


(27) 


From  figure  4  (E)  and  (F),  and  the  fact  that  the  Cl  at  1 ,  2,  4,  and  6 
are  simply  related  by  rotations,  it  follows  that 


I  2*q(  j  )SjC10(Sj  )Cn+2^m+q(Rj  )  -  4Zi<2iq(  1  )cn  +  2,m+q(  1  )<5m,4p  ' 

j 

where  the  prime  on  the  sum  means  j  =  3  and  5  are  missing.  For 
have  the  relation, 


O*  (5 )C  „  (5)  =  (-!)n+1pf  ( 3 )C  .  (3) 

q  n  +  2,m-Kj  1q  n+2,-m-q 


(28) 


5,  we 


and 


V3’  ‘  S.0(3,Sq.O 

from  the  symmetry  of  the  Cl-  site.  Thus, 

1  2*q(  j  'sjC1C)(sj  )Cn+2,m+q(Rj  )  =  2z2®*q  ^  3  ^Cn+2 , 0 (  3  '  ^mO  ^qO^n ,  2k+1  ' 

j-3,5  t29) 

where  k  =  0,  1 ,  2  .  •  . 

Finally,  using  equations  (28)  and  (29)  in  equation  (27),  we 

obtain 
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A,0(n  =  K,t2  ^6  [OT/DC^d)  +  ^,,(1)C*1(1  )]  +  32-,Qin(3)Cin(3)!  , 


11  31 


11  31 


2  10  30 


A30(,)  =  K3<4  /15  l^i(,)c5i(,)  +  ^il(1)C5i(,)l  +  ,0z2giO<3)So(3)f  ' 
A5q(1)  =  K5{12  /7  [P*,(  1  )C71  (  1  )  -  Qnd)C*l(1)'  +  21z2Ql0(3)C70(3)l  , 


A54(1)  =  K5{2  /3°  [Vll  Q*1(1)C75(1)  -  Q1i(1)c73(1)]}  ,  (30) 

A_(1)  =  K?{24  /5  iQ^mc^M)  +  QntDC^d)]  +  36z  ^  Q<  3  )CgQ(  3 ) }  , 
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A740  )  =  K7l4[/273  ^,d)C95(1)  -  3/5  Q,  ,  (  1  )Cg3  ( 1  )  ]  }  , 


where 


2  „  ,  n+3 

-e  4zi/R 


Also,  because  of  the  symmetry  of  the  Cl  site, 


Qui(1)  =  f  7=qio(3) 

v  2 


The  results  given  for  the  Anm  in  equations  (26)  and  (30)  are 
only  for  the  symmetry  coordinates  shown  in  figure  4  (E)  and  (F),  For 
the  true  normal  modes,  we  should  replace  z1  and  Zj  in  the  above 
by  z”  and  zj’,  where 

z"  =  z1  cos  i|/  +  z^  sin  , 

z'2  -  z  2  cos  i^tz^sin^  ,  (31) 


for  the  normal  modes  of  vibration  for  the  two  r.(t.  )  modes.  Two  values 

4  1  u 

of  ij;  in  equation  (31)  are  determined  from  the  force  constants  in  the 

vibration  Hamiltonian.  It  can  be  shown  that  the  effective  A  (k)  for 

the  two  normal  modes  of  the  r  (t  )  are  "  * 

4  1  u 

A"  (k)  =  A  (k)  cos  ill  +  A1  (k)  sin  1/ 
nm  nm  nm 

and 

A"  (k)  =  -A  (k)  sin  il  +  A’  (k)  cos  ih  ,  (32) 

nm  nm  nm 

where  the  A  (k)  are  given  by  equations  (26)  or  (30)  and  the  A'  (k)  are 
.  run  ,  .  ,  ,  _  ,  nm 

given  by  the  same  equations  but  with  z1  replaced  by  zj,  etc. 
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5.3 


Computation  of  A  (rt)  for  r  (t  )  Mode 
run  5  2u 


As  was  stated  earlier,  the  importance  of  the  vibronic  inter¬ 
action  is  to  mix  a  and  .1  states  with  the  f  states  to  allow  normally 
forbidden  electric-dipole  matrix  elements  of  those  ions  in  sites  with 
inversion  symmetry.  The  only  vibronics  that  can  contribute  to  the 
electric-dipole  transition  are  those  with  odd  parity:  r  1 .  The  refer¬ 
ence  displacement  relative  to  the  rare-earth  ion,  z1  in  equations  (21), 
(22),  (26),  and  (30)  is  simply  proportional  to  the  amplitude  of  the 

normal  mode,  say  Qk .  If  we  consider  transitions  where  a  phonon  is 
created,  then 

<l|Qk|0>=  (ti/2u!kMk)l/2  ,  (33) 


where  0Jk  is  the  frequency  of  the  vibronic  and  is  the  reduced  mass  for 
that  mode.  For  the  rg(t2u)  mode  (referred  to  as  Vg  elsewhere15),  the 
reduced  mass  for  the  (ErCl)-2  skeleton  is  simply  the  mass  of  the  chlo¬ 
rine,  and  equation  (33)  gives 


<1  I’3ki0>  =  0.073  A  ,  (34) 

where  we  have  taken*  the  energy  of  the  T^lt^)  mode  as  87  an’1 ,  as  given 
by  Hasan  and  Richardson.15  The  Cnm(j)  that  occur  in  the  expressions  for 
the  Anm(k)  can  be  obtained  from  Karayianis  and  Morrison,27  and  the 
values  of  R  can  be  obtained  from  the  x-ray  data.  Table  3  gives  the 
results  of  the  calculation  using  the  point-charge  contribution,  equation 
(21),  for  the  Anm ,  arising  from  the  r^(t2u)  mode  of  vibration.  It  is 
interesting  to  compare  the  results  given  in  table  3  to  the  odd-fold 
crystal-f ield  components  of  LiYF^,  which  does  not  have  an  inversion 
center.  For  Li YF4 ,  A32  =  657  -  667i ,  a52  =  -2671  -  59i ,  A72  =  7  +  1 4i , 

and  A7g  =  254  +  45i  (all  Anm  in  units  of  cm-1/An,  as  in  table  3).  The 

magnitude  of  the  Ann]  for  LiYF^  relative  to  the  effective  fields  due  to 
the  vibronic  interaction  in  CsjNaErClg  indicates,  at  least  for  the 
mode  r5(t2  ),  that  we  should  expect  the  latter  to  have  a  much  smaller 
line  strength.  Using  the  results  of  table  3,  we  calculated  the  intensi¬ 
ties  of  the  r3(t2u)  mode  for  the  monopole  interaction  of  Er^+  in 
Cs2NaErClg  for  all  the  Er^+  levels  below  27,000  cm-1.  The  results  of 
this  calculation  given  in  table  4  (p  28)  are  only  for  the  transition 
from  the  Tg  ground  state  to  the  various  upper  levels.  Thus,  the  calcu¬ 
lations  should  be  compared  to  experimental  data  taken  at  low  tempera¬ 
ture.  In  general,  for  all  the  multiplets,  the  intensities  of  the 

vibronics  are  approximately  equal  to  the  magnetic  intensities.  For 
the  ^Ii5//2  +  4*13/2  transition.  our  calculations  agree  qualitatively 

with  the  calculations  of  Hasan  and  Richardson.16  That  is,  the  vibronic 
intensity  is  smaller  than  the  magnetic-dipole  intensity  and  the 

Fg  ♦  Tg  magnetic-dipole  intensity  is  predicted  to  be  the  largest  (see 
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fig.  1,  Hasan  and  Richardson‘S).  The  calculation  for  the  ^1^5/2  *  ^F9/2 
tranisitions  indicates  that  all  the  vibronics  have  approximately  the 
same  intensity  as  the  corresponding  magnetic-dipole  transitions,  and 
this  corresponds  to  what  is  observed  experimentally  (fig.  3,  Hasan  and 
Richardson  :  8  ) .  a  detailed  comparison  of  calculated  and  experimental 
results  will  have  to  wait  until  the  dipole  contribution  and  perhaps  trie 
self-induced  interactions,  equation  (8),  are  included  in  the  effective 
odd-fold  field  components. 


TABLE  3.  P 1IMT- CHARGE  CONThliV'TI JN  To  OiJlJ  n  A  (cn"*/An); 

nn  ,  v. 

EFFECTIVE  CRYSTAL  -I ELDS  EOS  (  t,  )  MODE  JE  Cs^NaErClg  " 
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6.  CONCLUSION 

The  results  of  the  simple  calculation  of  the  intensity  of  the 
vibronics  of  Er3+  in  CSjNaErCl^  is  rather  encouraging  in  that  the  quali¬ 
tative  agreement  with  experiment  appears  good.  Certainly  the  calcula¬ 
tions  indicate  that  a  further  refinement  in  the  theory  of  the  vibration¬ 
al  interaction  and  the  inclusion  of  an  intensity  calculation  for  the  two 
t1u  vibrational  modes  is  called  for.  The  calculations  here  only  con¬ 
sidered  the  vibrations  of  the  (ErCi^)-3  complex  and  therefore  are  only 
approximate;  further  refinements  should  also  be  carried  out  with  the  use 
of  the  results  given  by  Lentz18  for  all  of  the  modes  of  the  unit  cell. 
Such  a  calculation  could  easily  be  carried  out  using  the  results  of 
equation  (17)  for  the  vibrational-electronic  interaction.  A  con¬ 
siderable  amount  of  experimental  data  exist  for  comparison  with  the 
calculation  for  many  R"*  +  in  CsjNaRCl^.  The  resulting  theory  can  then  be 
applied  to  the  possibility  of  constructing  a  tunable  laser  using  the 
isomorphic  material  Rb2NaRFg»  which  is  more  stable.  The  theory  is  also 
applicable  to  the  development  of  high-frequency  phonon  generation  (using 
the  lattice  modes  of  Lentz18).  Further,  the  theory  is  applicable  to  the 
generation  of  submillimeter  electromagnetic  waves  by  other  means  such  as 
polaritons  and  tuning  by  an  external  magnetic  field.  By  a  slight  exten¬ 
sion  of  the  theory,  the  possibility  of  building  a  tunable  green  laser 
using  Ce^+  in  RbjNaYF^  can  be  investigated. 
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TABLE  4 
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Kn^r-jy 

1  Energy*3 

Energy1' 

..nl 
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Phys.  7J_  (1980),  2580  (ref  25);  the  centroids  have  not  been  adjusted  for  best 
f  it . 

^Experimental ,  Hasan  and  Richer  Ison,  MoLec .  Phys.  45_  (1982),  1299  (re:  In); 
only  IR  labels  that  differ  from  a  are  listed.  The  levels  in  brackets  were  not 
observed  but  were  inferred  by  calculation  done  in  ref  1(1. 

cExperimental ,  Jezowska-Tr  zebiatowska  et  al,  Chen.  Phys.  50_  (1980),  209 

(ref  15);  IR  are  labeled  when  they  disagree  with  a. 


The  numbers  in  parenthesis  in  the  last  two  columns  are  powers  of  ten. 
IR--i rreducible  representation  of  the  cubic  group. 

,5mi-  -magnetic-dipole  strength. 


.sZ1-  -electric-dipole  strength  for 
u6 


v,  line. 
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